Discretely Exact Derivatives for Hyperbolic PDE-Constrained Optimization Problems Discretized by the Discontinuous Galerkin Method

نویسندگان

  • Lucas C. Wilcox
  • Georg Stadler
  • Tan Bui-Thanh
  • Omar Ghattas
چکیده

This paper discusses the computation of derivatives for optimization problems governed by linear hyperbolic systems of partial differential equations (PDEs) that are discretized by the discontinuousGalerkin (dG)method.An efficient and accurate computation of these derivatives is important, for instance, in inverse problems and optimal control problems. This computation is usually based on an adjoint PDE system, and the question addressed in this paper is how the discretization of this adjoint system should relate to the dG discretization of the hyperbolic state equation. Adjoint-based derivatives can either be computed before or after discretization; these two options are often referred to as the optimize-then-discretize and discretize-then-optimize approaches. We discuss the relation between these two options for dG discretizations in space and Runge–Kutta time integration. The influence of different dG formulations and of numerical quadrature is discussed. Discretely exact discretizations for several hyperbolic optimization problems are derived, including the advection equation, Maxwell’s equations and the coupled elastic-acoustic wave equation. We find that the discrete adjoint equation inherits a natural dG discretization from the discretization of the state equation and that the expressions for the discretely exact gradient often have to take into This document has been approved for public release; its distribution is unlimited. Lucas C. Wilcox (B) Department of Applied Mathematics, Naval Postgraduate School, Monterey, CA, USA e-mail: [email protected] G. Stadler · T. Bui-Thanh · O. Ghattas Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX, USA T. Bui-Thanh Department of Aerospace Engineering and Engineering Mechanics, The University of Texas at Austin, Austin, TX, USA O. Ghattas Departments of Mechanical Engineering and Jackson School of Geosciences, The University of Texas at Austin, Austin, TX, USA

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Investigation of Commutative Properties of Discontinuous Galerkin Methods in PDE Constrained Optimal Control Problems

The aim of this paper is to investigate commutative properties of a large family of discontinuous Galerkin (DG) methods applied to optimal control problems governed by the advection-diffusion equations. To compute numerical solutions of PDE constrained optimal control problems there are two main approaches: optimize-then-discretize and discretize-then-optimize. These two approaches do not alway...

متن کامل

A Posteriori Error Estimation for Discontinuous Galerkin Approximations of Hyperbolic Systems

This artic!e considers a posteriori error estimation of specified functionals for first-order systems of conservation laws discretized using the discontinuous Galerkin (DG) finite element method. Using duality techniques, we derive exact error representation formulas for both linear and nonlinear functionals given an associated bilinear or nonlinear variational form. _Veighted residual approxim...

متن کامل

A Posteriori Error Estimation for Discontinuos Galerkin Approximations of Hyperbolic Systems

This article considers a posteriori error estimation of specified functionals for first-order systems of conservation laws discretized using the discontinuous Galerkin (DG) finite element method. Using duality techniques, we derive exact error representation formulas for both linear and nonlinear functionals given an associated bilinear or nonlinear variational form. Weighted residual approxima...

متن کامل

Quadrature-free Implementation of the Discontinuous Galerkin Method for Hyperbolic Equations

A discontinuous Galerkin formulation that avoids the use of discrete quadrature formulas is described and applied to linear and nonlinear test problems in one and two space dimensions. This approach requires less computational time and storage than conventional implementations but preserves the compactness and robustness inherent to the discontinuous Galerkin method. Test problems include both ...

متن کامل

Discontinuous Galerkin methods for first-order hyperbolic problems

In this paper we consider discontinuous Galerkin (DG) finite element approximations of a model scalar linear hyperbolic equation. We show that in order to ensure continuous stabilization of the method it suffices to add a jump-penalty-term to the discretized equation. In particular, the method does not require upwinding in the usual sense. For a specific value of the penalty parameter we recove...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • J. Sci. Comput.

دوره 63  شماره 

صفحات  -

تاریخ انتشار 2015